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DUAL GRAPHS AND MODIFIED BARLOW-BASS RESISTANCE ESTIMATES FOR 
REPEATED BARYCENTRIC SUBDIVISIONS 

ANTONI BRZOSKA, DANIEL J. KELLEHER, HUGO PANZO, AND ALEXANDER TEPLYAEV 


Abstract. We prove Barlow-Bass type resistance estimates for two random walks associated with 
repeated barycentric subdivisions of a triangle. If the random walk jumps between the centers of 
triangles in the subdivision that have common sides, the resistance scales as a power of a constant 
p which is theoretically estimated to be in the interval 5/4 ^ p ^ 3/2, with a numerical estimate 
p « 1.306. This corresponds to the theoretical estimate of spectral dimension ds between 1.63 
and 1.77, with a numerical estimate ds ~ 1.74. On the other hand, if the random walk jumps 
between the corners of triangles in the subdivision, then the resistance scales as a power of a constant 
p T = 1 /p, which is theoretically estimated to be in the interval 2/3 ^ p T ^ 4/5. This corresponds 
to the spectral dimension between 2.28 and 2.38. The Barlow-Bass method for proving the resistance 
scaling for Sierpinski Carpets has to be modified considerably to achieve the results in the current 
work. 
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1. Introduction 


The repeated barycentric subdivision of a simplex is a classical and fundamental notion that ap¬ 
pears frequently in algebraic topology, see lH8l and references therein]. Recently it was considered 
from probabilistic point of view in [|9i [12l HQ, TO [361 . Understanding how resistance scales on 
finite approximating graphs is the first step to developing analysis on fractals and fractal-like struc¬ 
tures, such as self-similar graphs and groups, see [0 [1511221 :34, 23 and references therein]. For 
finitely ramified post-critically finite fractals, including nested fractals, the resistance scales by the 
same factor between any two levels of approximating graphs. Kigami, Linsdtrpm et al [24, 25, 30] 
used this fact to prove the existence of a Dirichlet form on the limiting fractal structures. In the 
infinitely ramified case, resistance estimates are more difficult to obtain, but are just as important 
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Figure i. Barycentric subdivision of a 2-simplex, the graphs Gj n G[ and G\. 

to understanding diffusions on fractals. Barlow and Bass 02 [2, 3HU proved such estimates for the 
Sierpinski carpet and its generalizations. These techniques were extended to understanding resis¬ 
tance estimates between more complicated regions of the Sierpinski carpet, see lf32l . The paper 
Il26l provides another technique for proving the existence of Dirichlet forms on non-finitely rami¬ 
fied self-similar fractals, which estimates the parameter p by studying the Poincare inequalities on 
the the approximating graphs of the fractals. 

Our work further develops existing techniques to obtain resistance scaling estimates for the 1- 
skeleton of n-timcs iterated barycentric subdivisions of a triangle which we will denote G^, and 
its weak dual the hexacarpet (introduced in | J6| ), which we will denote G„ . Note that in the current 
work, hexacarpet graphs will mostly be used to refer to G%, which is a modification of G n by adding 
a set of “boundary” verteces. 

If RT and R n is the resistance between the appropriate boundaries in Gfi and Gfi respectively 
(see Figures [l] [2[ [3} |4]), then we prove that the resistance Rfi and R n scale by constants p T and p 
respectively, and obtain estimates on these constants. Our main result is the following theorem. 

Theorem 1.1. The resistances across graphs Gj r and G% (defined in Section [2]) are reciprocals, 
that is 

Rn = 1 / R n 

and the asymptotic limits 

log p 1 = lim — log RIL and log p = lim — log R n 

n—>00 Tl n—>00 n 

exist. Further, p T = 1/p and satisfies the estimates 

2/3 < p T < 4/5 and 5/4 <p< 3/2. 

These estimates agree with the numerical experiments from [j6l, which suggest that there exists 
a limiting Dirichlet form on these fractals and estimates p « 1.306, and hence p T ~ 0.7655. 

Dirichlet forms have many applications in geometry, analysis and probability, see [0 [8] H4J. In 
general terms, a Dirichlet form on a fractal is a bilinear form which is analogous to the classic 
Dirichlet energy on R d given by 

<?(/) = f IV/I 2 dx. 

A Dirichlet form is equivalent, in a certain sense, to a symmetric Markov process on the base 
space. The potential theoretic properties of the Dirichlet form have implications for this stochastic 
process. In particular, the resistance between two boundary sets is related to to the crossing times 
of these sets. In our case, on the 1-skeleton G^, the Markov process jumps between corners of 
the triangles in the subdivision. Our theoretical estimates correspond to a process with limiting 
spectral dimension between 2.28 and 2.38. The Markov process on the hexacarpet graphs (which 
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Figure 2. Adjacency (dual) graph G 2 and the barycentric subdivision graph C' 2 
of a 2-simplex. 

are denoted by G and G H later on) corresponds to a random walk which jumps between the centers 
of these triangles, with spectral dimension between 1.63 and 1.77 (« 1.74 using the numerical 
estimates in 0). This is a substantial difference implying, in particular, that one Markov process 
is not recurrent, while the other is recurrent. From the point of view of fractal analysis, our results 
suggest that the corresponding self-similar diffusion is not unique, unlike IfTTl l4l. The long term 
motivation for our work comes from 120, 19, 2H [271 and related papers, and especially from the 
works on the heat kernel estimates [[Hill 4.HUES, 28.1291 

In the discrete setting, the Dirichlet form is the graph energy. In this case the resistance between 
two sets is determined using KirchhofFs laws. In the appendix in Section[5]we collect basic results 
about graph energies and resistance which we need throughout the paper. For a more thorough 
introduction to these topics, one can see, for example, llT3l l3~ll. 

This paper is organized as follows. Section [2] lays out the formal definitions and shows how to 
take advantage of the duality of the graphs G' 77 and G[ to prove that R n = 1/R[. In Section [ 3 ] we 
first prove sub-multiplicative estimates R m+n < cR m R n for some constant c independent of m and 
n. This is done in a fashion generalized from 0]|. Fekete’s theorem implies that the limits p and p T 
exist. To prove that these limits are finite, in section [4] we establish upper and lower estimates on 
Rj t and R n . This is done by comparing G T n and G\[ to subgraphs and quotient graphs respectively. 

2. Barycentric subdivision, the hexacarpet and the resistance problem 

We define the 2-simplicial complexes T n = (E°, E',. E’;J where E( are the /'-simplexes of the 
complex, starting with a 2-simplex (a triangle) T 0 with O-simplexes (vertexes) Efj = {po,Pi ■ P 2 }, 
1-simplexes (edges) E.J = {[puPj]}^, and 2-simplex (triangle) E 5 = {\po,pi,p 2 ]}. Here, if 
q 0 , (p, q 2 e E°, then \q {) , <p will refer to the 1-simplex with q 0 and q 1 as endpoints (which may 
or may not be in E 7 ), and similarly for [g 0 , <p , q 2 \. We will only be considering simple simplicial 
complexes (without multiple edges/triangles) and we shall not be considering orientation. Thus 
|i/o- <7i] = \ ( l\■ 7 o| determines a unique 1-simplex. We shall also use the notation < to denote con¬ 
tainment of simplexes, i.e. q 0 < [qo, qi] < [ 90 ) 91 ) 92 ]- 
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Figure 3 . Adjacency (dual) graphs G 4 and Gj for barycentric subdivision of a 2-simplex. 

T n+ i is defined inductively from T„ by barycentric subdivision, pictured in figure [I] That is 
E° +1 is E° along with the barycenters of simplexes in E' ( , i = 1,2, which we shall refer to as 
b(e) for e G E' n , i = 1, 2. 1- and 2-simplexes of T n+] are formed from barycenters of nested 
simplexes. More concretely put, elements of E* +1 are either of the form [g 0 , b([g 0 , q } ])] where 
Qo,Qi e E° with [q 0 , qi \ G E*, [q 0 , b([g 0 , qi, g 2 ])], or [b([g 0 , 9i]), KWp^vS)} where q 0 , qi , q 2 G 
E° with [g 0 , gi, g 2 ] e El, and elements of E* are of the form [g 0 , b([g 0 , gi]), b([g 0 , q u g 2 ])] where 
qo,qi,Q 2 e E° with [go, gi] G E, 1 , and [g 0 , gi, g 2 ] G E;. 

Definition 2.1. We define the graph G 4l = (V^, Ej t ) with vertex set Vj = E° and edge relation 

q ~ T q' if [g, q'\ G E' ( . We will refer to this as the 1-skeleton of the T n . 

Definition 2.2. Following |0, we define the graph G n = (V n ,E n ) with vertex set V n = E 2 n and 
edge relation [g 0 , gi, g] ~ [go, gi, g'|. That is, the vertexes of G„ are the 2-simplexes of T n and they 
are connected by an edge if these simplexes share a 1-simplex. 

Remark 1. Gj L is classically known to be a planar graph, as seen in figure [lj although throughout 
this work we will refer to the hexagonal embedding from figure [4] more often. With either of these 
embeddings, G n is the weak planar dual, that is each of its vertexes correspond to a plane region 
carved out by the embedding of G ' n with the exception of the unbounded component. 

We will need explicit names for the elements of E^ = {[> 0 , Pi , P 2 , Po , Pi , P 2 , I> } where {p 0 , pi , p 2 } = 
E[j as above, p\ = b ([pi,Pi']), i' = i + 1 mod 3, andp' = b(\p 0 ,pi,p 2 }). 

This is convenient for recursively defining functions on T n . Define self-similarity maps F t : 

T n -A T n+ i for i = 0,1,..., 5, which are defined on T 0 by Fi(p 0 ) = p', F^p 1 ) = pp/ 2 ], and 

Ffa) = P[i/ 2 \ > where the index is taken mod 3. F, is extended to T n by the relations Fj([g 0 , gi]) = 
[Fi(g 0 ),Fi(gi)], Fi([g 0 ,gi,g 2 ]) = [Fj(g 0 ), F^), Fj(g 2 )] and b o F t = F { o b. If w = w 4 w 2 ■ ■ ■ w k 
is a word in (0,1,..., 5} fc , then we define F w : T^ —> Tj l+k by F w := F Wl o F W2 o ■ ■ ■ o F Wk . 

Fi as a function from G^ Gf l+1 or G n —> G n+ 1 is a graph homomorphism, and G„ +1 = 

U f =0 Fi(Gn). However, this is not true for G n and G n+1 , since not every edge of G n+ 1 is covered 
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Figure 4 . A^)* and B^ H on the hexagonal embedding of G T J H . 

by Fi(G n ) for some i = 0,1,..., 5. We want to take advantage of this self-similarity throughout 
the current work, thus we shall define a modified hexacarpet graph. 

Definition 2.3. The (modified) hexacarpet graph G// is defined to have vertex set 

V H := E 2 U E 1 

n n n 

where adjacency is determined by [qo,qi,q\ ~h [?0)9i]> he., e G E/ is connected to / G E 2 if 

e < f ■ 

We define the graph energy and energy dissipation on these graphs as in the appendix in section 
[ 5 } so it is sufficient to define the conductance of edges. For Gj x , 

f 0 if [q, q'] i E l n 

c Iq' = \ 1 if there exists ^0 7 ^ Qo with [q, q', q 0 ], [q, q', q 1 ^ G E 2 
I 2 if there exist only one q 0 such that \q, q', q 0 ] G E 2 . 

We take to be the graph energy defined with the above conductance. The advantage of these 
conductance values is the resulting self-similarity relation 

-Ci(/) = E'C(/°a). 

i=0 

Similarly if we define c^ q , = 1/2 if q q' and 0 otherwise, then the resulting energy function 
satisfies the following relation 

i=0 

Both of these relations are also true for energy dissipation of functions on edges of these graphs. 
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It will often be useful to think of these graphs as embedded in the plane [R 2 . It is typical to think 
of T n as a subdivided triangle in the plane, but we shall embed it as a hexagon, as T n , n > 0, 
has symmetry group D 6 , the dihedral group on 6 elements. As such, for n > 0 we define a map 

F T :E^ = by F T (p') = (0,0), 

F T { pk ) = (cos(2/c7t/3), sin(2A;7r/3)) and F J (p' k ) = (cos((2fc + 1)7t/6), sin((2fc + l)7r/6)) 


for k — 0,1, 2. Thus is mapped to the corners and midpoint of a regular hexagon centered at 
(0, 0), see figure |4j We extend to E° by taking averages: 


F T ob([q 0 , qi }) 


F T (q 0 )+ F T (g i) 

2 


and F r ob([q 0 ,q 1 ,q 2 \) 


F T (q 0 ) + F T ( qi ) + F T (q 2 ) 
3 


We embed by the map F u : Vff —» R 2 by F H = F T o b. Thus the vertexes of the embedded 
G\[ are the centers of the triangles and edges of the embedding of Gj t . 

If e = [q 0 , q\\ G E/ then we define the geometric realization of e , |e|, to be the convex hull of 
F T (p 0 ) and F T (p 1 ). That is 


|e| = {6 0 F T (p 0 ) + d 1 F T (p 1 ) : 9 0 + 6 1 = l} . 

Similarly, if / = [q 0 , q \, q 2 ] G E 2 , then |/| is defined to be the convex hull of F T (p 0 ), F T (pf) and 

F T (p 2 ). 

To define the resistance problem on these graphs, we shall need to define the boundary of these 
graphs. N c |/|, for e G E l k and / G E' ( for k < n, i < j if the geometric realization of e is a 
subset of the geometric realization of/, e.g. |b(e)| C |e| or \[q 0 , b([g 0 , qi])]\ C [go, q \, 92 ] |- 
Define L[ n,T \ resp. L[ n,H \ i = 0, 1,... 5 to be the vertexes q G V r f, resp. , such that 
|g| C | [pi/ 2 , p\/ 2 ) | if i is even, and the set of \q\ C | \pj,p'u-iy 2 \ \ where j = (i + l)/2 mod 3 if i is 
odd. We will suppress arguments of the superscript when there is no danger of confusion. 


Definition 2.4. Define U L ( ] n ’ H) and = lJ : !'' ir> U L^ l ' H> . Further, define R n to 

be the effective resistance with respect to <§k^ between A^ n) and . 


Definition 2.5. Define Aj^ = il^' 1 and B, r nj = L ( F T> . Further, define /?/ to be the effective 
resistance with respect to S. [/ between Aj n) and BF. 


Theorem 2.1. The resistances IT,, and /i'/ are related by 


< = T5~- 

-Ll'n, 


Proof. The main tool of thi s pro of is proposition 9.4 from PTTl . Define the weighted graph Gj t = 
(Vj, Ef ) as in Proposition 


5.3 


where Vj is V r ; 


modulo the relation which identifies all elements 


of y\J n) and Bj n ^ into two vertexes called of and b T respectively. From proposition |5.3| effective 
resistance between a T and b T with respect to Aj is /?./. 

Further we define G r J by, not only identifying Af . and BF into single vertexes a H and b H 


respectively, but also to replace the sequential edges of the form [q 0l q \, q] [q 0 , qf [q 0l g,, q'\ 
using Kirchoff’s laws. Thus the sequential connections with resistance 1/2 are replaced with one 
connection [q 0 , q u q] [q 0 , qi,q'] with resistance 1. It is easy to see that the resistance between 
a H and b H is R n . Also, remove the vertexes contained in Aj n) and Bj n ^ and associated edges — 
since these vertexes are connected to the graph by only 1 edge,removing them has no impact on the 
resistance. 
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Figure 5. G 2 and G ' 2 without the additional edges. 

If we define the graphs (G^V an d (GnY t0 the G% and G' n with an additional edge connecting 
a H to b H and a 1 to b 1 respectively, then (G^) ' 1 is the planar dual of (Gj % ) t (see figure [ 5 ]), and thus 
by Proposition 9.4 in Oil . R n = 1 /R^. □ 


3 . Estimates 

This section is dedicated to the proof of 

Theorem 3.1. R m+n < | R n R m , and equivalently R^ n+n > \R^ n R^ l . 

This theorem implies that there are constants c 0 and ci such that c 0 + log R[ and ci + log R n 
are superadditive/subadditive positive sequences. Thus the limits exist, and we have the following 
corollary. 

Corollary 3.1. The limits 

lim — log R'E = log p T and lim — log R n = log p 

n—>00 Tl n—>• 00 Tl 

exist. 

Note that this corollary does not rule out the possibility that p = 0 or 00. In section |4j we 
establish positive upper and lower estimates on p and p T . 

We give two proofs of theorem |3Tj the constants differ in the proofs, but the exact value of the 
constant does not effect the existence of p and p T . This establishes the result independently of the 
duality of the graphs — indeed variations on these proof work for different choices of boundary. 

One version proves the upper estimate on R m+n directly, and uses flows on G r J . The direct proof 
of the lower estimate on Rj n+n is proven using potentials on Gj t . The two proofs mirror the upper 
and lower bounds for the resistance of the pre-carpet approximations for the Sierpinski carpet in 
m. The two versions of our proof highlight the importance of duality in proving Barlow-Bass style 
resistance estimates, and suggest possible directions in which their proofs can be generalized. 
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t 



Figure 6 . The transformation from the flow I n (left) to the flow Hq 2 (right) 


3.1. Upper estimate and flows on G^. Consider the flow I n which is the minimizing flow on 
from A?\ to . The general idea is to construct a flow on G n+m which has global structure 
resembling I m but on F^G^ for u; € {0,1,..., 5} m has structure which is built from I n . 

Recall that A^ n) = l} ( " } U L 1 ]"' 1 , then I n is a flow from /\^ n) to define // ( " 2 = I 71 when re¬ 
stricted to the half of the G r J which connects vertexes contained in [p 0 ,p' 0 ,p'}, [p' 0 , pi, p '\, [p i. p\ ■ p'] ■ 
On the other half, Hq 2 is I n composed with the symmetry which exchanges po with p \, p' 0 with pi, 
and p 2 with p' 2 and is extended to the rest of by convex combinations. The construction of Hq 2 
is depicted in figure [6| Using the embedding F H . this symmetry is the reflection through the line 
which makes a tt/2 angle with the x-axis. Because I n subjected to a 7r-rotation is — I n , Hq 2 is a 
flow. 

Thus Hq 2 is a flow on G% between U Lj T,) and U L^\ with flux 1. is the flow 
from lJ ( ['U l|"' ) and lJ 2 ] U with flux 1 obtained from Hq 2 by applying the appropriate 
symmetry. Since the values of (resp. Hq 2 ) are just a rearrangement of those in I n , then 

E(H? n ) = E{Hq 2 ) = E{I r *) = R n for all i ^ j. 

We shall consider as the set of Y-networks, centered at x 6 C V^j. Takea 0 (x), ai(x), a 2 (x) 

to refer to the outward flow of I m restricted to each of these edges in the Y-network associated to x 
with orientation such that a 0 (x) is of a different sign then ai(x) and a 2 (x). i.e. if a 0 (x) < 0, then 
a\(x), a 2 (x) >, if a 0 (x) > 0, ai(x), a 2 (x) < 0, and in the case when a 0 (x) = a i(x) = a 2 (x) = 0 
then the choice is arbitrary. Notice that 

E(I m ) = ^ ^ (a 0 (x) 2 + ai(x) 2 + a 2 (x) 2 ) = R m 

Take F U G" C G^ +m such that FG\j> {) . p\. p 2 ]) = x to be the subgraph of G^ +m isomorphic 
to G r J corresponding to x G E 2 ra . We shall assume that the labeling of the sides FT( L 2 U L^) be 
contained in the edge which corresponds to the values a\ (x), and F U (L2 ULg) corresponds to a 2 (x). 
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Figure 7 . The function u, v and w. 

The set of subgraphs {RuiGn} xeE2 wg | 0 x ~y n . cover G^ +n and define the flow J on G ^ +n by its 
values on these subgraphs 

JoF u = a 1 (x)HS 1 + a 2 (x)HS 2 . 

J is a well defined flow because Hq 2 is obtained by a reflection of which I n is symmetric with 
respect to, so aHfi n (y, z) + bHRfiy, z) = (a + b)I n (y, z) for all y e Lq U Lfi. 

Now we see, 

E(J)= £ E(a 1 (x)HS l +a 2 (x)HS l ) 

= ]T (a 1 (x) 2 E(H^) + a 2 (x) 2 E(H^) + 2a 1 (x)a 2 (x)E(H^,H^)) 

x £ E ?n 

< £ + a ^) 2 E (ff„” 2 ) + 2a 1 (x)a,(x)E(HS 1 ) 1/2 E(HS 2 ) 1/2 ) 

= (ai(x) + a 2 (x)) 2 E(4) = R n ^ a 0 (x) 2 

xSE^ xeE2 n 

2 x 4 

< Rn- 2^ (ao(^) 2 + ai(x) 2 + a 2 (x) 2 ) = -R m R n . 

Note that the first inequality holds because of Cauchy-Schwartz inequality and the fact that, by 
our labeling convention, ai(x)a 2 (x) > 0, and the last inequality holds because a 2 (x) + al(x) > 
(ai(x) + a 2 (x)) 2 /2 = ao(x)/2. 

3.2. Lower estimate and potentials on Gj t . Let o n be the harmonic potential on G' n with bound¬ 
ary values 0 on AR and 1 on BR . On G^_ 1 , define u — (f> n o F 0 , v — (fi n o Fi and w = fi n o F 5 . 
<fi n o F 2 , (fi n o F 3 , and <fi n o F 4 can be written in terms of u, v, w and the constant 1 function as 
illustrated in figure |T| 
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Notice that the function w is v o a where a : Gf_ , —» Gf_ , is the symmetry which exchanges 
Pq and p 2 , fixing p 0 , and is extended to the rest of Gf_ , by averages ( with respect to F T , a is the 
flip about the horizontal axis). Also, 


K l = tfitn) = 2(Cl («) + <Cl(«) + ClH) = 2<Cl(«) + 

Lemma 3.2. fu, n — w) = 0 

Proof. On one hand, Sff_ ff) = o cr) for all / because a is a graph isometry with respect 

to the conductances. However, the function u symmetric about the horizontal axis, i.e. u(x) = 
uoa(x), and v and w is anti-symmetric about this axis, i.e. (voa(x) — woa(x)) = (w(x) — v(x)). 
Thus 


<C-1 (u, V-w) = O a,v O (7 — w o a) = w - v) = V-w). 


□ 


For each x e E^_ 1 , i' = i + 1 mod 6 define a x = b(x) to be the barycenter of x. Also, define 
af,i = 0,1,..., 5 to be vertexes and barycenters of edges contained in x ordered in such at way 
that [of, of,] is an edge in E*, and that the vectors 


(v o F u (a x ),vo F u (a x ),v o F w (of,)) = (1/2, 0,1/2), 

{w o fG(a x ), w o F u (af),w o F u (a x ,)) = (1/2,1/2, 0), and 
(u o F u (a x ),u o F u (a*),uo F u (a x ,)) = (1/2, 0,0), 

forte G {0,1,..., F] n ^ 1 such that F u ,(\p 0 ,p l ,p 2 \) = x. Notice that for all x e E^_ l5 a x and a x are 
contained in E°, and that, for any function / on Gf. then 

€ U )= E E ((/ m -/(« r )) 2 -5(/«)-/(“?)) 2 


where i' = i + 1 mod 6. We now define a function f n+m on Gj l+rn such that / n+m | G T = f> m as 
follows: if F u : Gf_ , — y Gf n+n is the contraction mapping which takes Gf_ , to the \a x , of, off, 
then 


fn+m O F u = (2(j) m (a X ) - 4> m (of ) - 4> m (of,) )«■+ (0 m (of ) - «/) ) (v ~ w) + - (0 m (of ) + (j) m (of,) ) 1. 
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So 


&n+m{fn+m)= E C-l(/n+iB°^i) 
we{o,...,5} m+1 
5 

= E E^™^ - <^«) - ^(a?)) 2 ^-!^) 


mSE 


(2) i=0 


+ (&»(<»?) - i«))’C,(t> - m 


< ^(2^,(0') - A»(of) - ^(o?))X-i(«) 


xeE 


( 2 ) 1=0 


+ 4(0 m «O-0 m «))X T _ 1 (^ 


< 


E ]T 2((^(a*) - <M<,)) 2 + (<UO - 0 m «O) 2 )<Ci(«) 


xeE 


( 2 ) 1=0 


+ 4(0 m «,)-0 m «))X T _iW 

< 2 (2^_ 1 (u) + 4(f^_ 1 (w)) 


^ <=o V 7 




Thus i? m+n > R m R n /2. 


4. Estimating resistance using open and short circuits 
In this section we prove estimates on the constants p and p T to prove that they are finite. 

4.1 . Upper bound on p by removing edges. 

Proposition 4.1. R n < (3/2)”, and thus p < 3/2. 

Proof. Figure [8] shows G\ = (V// Ef), which is obtained from Gf by removing all edges which 
have the vertexes contained in |[pg,p']| and [//. p' \ . Figure [s] also illustrates how six G n graphs 
can be glued together to form one G n+ 1 graph. By induction, we see that each G n is made up of 2 n 
paths between L U L ( -f’ H ' > and U L^’ H \ 

The lengths of all the paths from Lq 1,h '* U L^ l,H ' > and U L^ 1 '^ in G n can be encoded in 

a sequence of 2 n integers. Let’s call this sequence l n and write l n = {l n ,i,ln, 2 , where 

l nj is the length of the path which has initial (or terminal) point j th closest to p 0 (in graph or 
Euclidean distance with our embedding). Since G n+ 1 is 6 copies of G n glued together, }_^ k=l l n ,k = 
6 SLi l n-i,k = 6 n , because h = 2,4. 
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The corresponding resistance R n between Lq L ' h ^ U L ( f n> and I3G' H) U L ( f Jr> of the G n graph 
is the resistance of 2 n paths connected in parallel so, by Kirchhoff’s laws, 

( 2 n \ ~ 

££) ■ 

Recalling that a; i—>• — is convex on (0, oo), this allows us to use Jensen’s inequality to write 



To obtain an upper bound on R n , we glue together 6 copies of G n -\ to produce a graph which 
has an edge set contained in Et* consisting of 2 n paths from A r J to 13 f. Then, using proposition 


to with 2 parallel connections of three sequential wires of resistance R n . Thus we have that 

Rn < ZRn-l/2 = (3/2)-. □ 

This also implies a lower bound on p T > 2/3. However, this could also have been achieved by 
considering the graph G' t) with vertex set Vf where x. y <6 Vj (the vertex set of G ' n ) are identified 
as one vertex if an edge connecting x and y was deleted in the construction of G n . One attains a 
graph as in figure [9], with “end” points P and 0, which correspond to the points in A T and I3 r . 

4.2. Lower bound on p by shorting graph. In this section we define a graph G%, such that V n 
is modulo an equivalence relation. Thus, resistance in G^ between two sets is less than the 
resistance between the fibers of these sets. 

Proposition 4.2. There is a constant c independent of n such that R n > c(5/4) n . Thus p > 5/4. 


5.2 and Kirchhoff’s laws in addition to the above argument, we obtain that we are connecting Af 
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Figure 9. 



Figure 10. Graph G 2 with short circuits. 


Proof. We define two vertexes in Vj 1 to be equivalent if they are both in contained in F L0 (\jj l ,p J ]) 
where j = i+ 1 mod 3 for some i and some u e {0,1,..., 5} fc . Figure[lo| shows G 2 ■ These graphs 
appeared in 051 . as an examples of non-p.c.f. Sierpinski gaskets, where it was determined that their 
resistance scaling factor is |. This implies that R n +\ = | R n and subsequently R n = Ri (§)” _1 . 
From this, and a gluing argument as in the previous subsection, it follows that R n > c(|) n for some 
c > 0 which is independent of n. □ 


Using a similar argument, G' n can be obtained by identifying points in graph approximations 
of another non-p.c.f. Sierpinski gasket which appears at the end of ||35|. and which is pictured in 
figure 11 The resistance between the corner points of these graphs is some constant times (4/5) n , 
and thus the resistance between the corner points of G' n is greater than this value. This proves that 
p T < 4/5 because including more points in the boundary decreases resistance, so the resistance 
between the corner points is greater than the resistance between A and B. 

Alternatively, if you connect the boundary points of the graph in figure [TO] to a A-network, it is 
dual to the network attained by connecting the corner points in the graph in figure[TT]to a Y -network. 
This explains why the resistances are reciprocal. 


5. Appendix: energy, potentials, and flows on graphs 

In this section we collect the basic definitions and results concerning graph energies and resis¬ 
tances which are used in the current work. For more detailed expositions on this subject see, for 
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instance, HD USED. Let G = (V, E) be a graph with finite vertex set V, and edge set E, which is 
a symmetric subset of V x V. We define the set £(V) = {/ : V —>■ R} to be the set of real valued 
functions on V which we will sometimes refer to as functions or potentials on the graph G. 

If G is a weighted graph there are conductances (weights) c pj] = c q . p for all p,q e V such that 
c P)q > 0 if (p, q) G E and c P)Q = 0 when (p, q) E. Resistances between points p and q for any 
(p, q) G E will be defined r m : = 1/c P)9 . To a graph with associated conductances, define the graph 
energy 

£ : £(V) x £(V) —>■ IR, such that &(f,g) = ^ c pM(p) ~ f(q))(g(p) ~ $(?)), 

(p,q)£E 

and write S’(f) := <£*(/, /). The vector space (:(V r ) is given the inner product 

</>0> = J2f{p)g(p)- 

p£V 

Functions (with orientation) on the edge set will be denoted 

£(E) = {J : V x V —> [R [ J(p, g) = —J{q,p), and J(p, q) — 0 if (p, q) ^ E} . 

Define the energy dissipation E : £{E) x £{E) —> IR to be the inner product 

E ( J , K ) = \ r p,q J (P^) K (P^) 

(p,q)£E 

on £{E). The discrete gradient V : £(V) —» £(E) is given by V/(p, g) = c P:q (f(p ) — /(g)), and 
discrete divergence div : £{E) —> £(V), defined by div J(p) = — J2 q ■ ( pq )eE •((/>■ <?)• Wc adopt 
the notation V and div because 

E(V/,J) = ^ ^ (f(p)-f(q))J(p,q) = '52 f(p) J (P^) = ~ (/ 5 divJ )> 

(p,q)£E p&V (p,q)£E 

and the weighted graph Laplacian A : £(V) -A ((E), defined by 

-A/(p) = c p>q {f{p ) - /(g)), 

(p,q)£E 

can be interpreted as A = div V. Further, it is elementary to check that the graph energy can be 
represented as 

£{f,g) = E{Vf,Vg)=-{g,Af). 

For A, B C V, J e ((E) is called a flow from A to B, if div J(p) = 0 for p ^ A U B. The flux 
of a flow J from A to B is defined by 

flux(A, B,J) = Y] div J{p) = ~Y^ div J{P)- 

P&A P eB 

These two expressions agree because, if J is a flow from A to B, 

div J(p) + ^ div J(p) = J^div J(p) = ^ ^ (J(p, q) + J{q,p)) = 0. 

p&A p&B pSV (Pi?)S-E 

A flow J from A to B is called a unit flow if flux (A, B, J ) = 1. Notice that if J is a flow from A 
to B, then 

E(V/, J) = — (/, div J) = — J2 f (p) div J(p) ■ 

p&aub 
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Figure 11 . Non-p.c.f Sierpinski gasket 

Additionally, for any flow J from /I to B where A and B are disjoint, and any function f\ A = 0 
and f \ B = 1, 

(/, div J) = f(P ) div J(p) = div J (P) = ~ fl ux(A J )• 

pSV p&B 

The effective resistance between sets A and B is defined by 

B)~ l = inf {S’(f) | /|a = 0, f\ B = 1} • 

Energy is minimized by the function 0 such that 0 | A = 0, 4>\b = 1 and A0(p) = 0 for p ^ A U B, 
and thus ^(0) = 1 / R(A. B ). We shall refer to such a function 0 as the harmonic function with 
boundary A and B. The only functions which satisfies A/ = 0 (i.e. harmonic without boundary) 
and f\ A uB = 0 is the constant 0 function, thus the 0 is unique. 

Proposition 5.1 (Thompson’s principle). There is a unique unit flow I which minimizes energy 
dissipation, and I = R(A, /i) Vg, where 0 is the harmonic function with boundary conditions 
0 \ A = 0 and 0 |b = 1. That is 

E (/) = inf {E( J) | J is a unit flow from AtoB}. 

Proof Because A0(p) = 0 for p 0 A U B, it is clear that V0 is a flow from A to B. We check 
flux(A, B, V0) = — div V0 = — E A 0(p) 

pen peB 

= - J] M^P) = - = ^(0) = S)" 1 . 

peB peu 

The third to last equality is because for p e V \ B, either 0(p) = 0 or A0(p) = 0. Thus / = 
f?(A, B) V0 is a unit flow. 

Further, if J was another unit flux flow, flux(A, B,J — I) = flux(A, B , J ) — flux(A, B , /) = 0, 
thus J = / + A' where A' is a 0 flux flow from A to B. Now, it is elementary to check 

E( J + K) = E(J) + E(A') + 2 R(A, B) E(V0, K) 

= E(J) + E(AT) + 2 R(A, B ) flux(A, B, K) 

= E(7) + E(AT) > E(J). 
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Thus / is the minimizing unit flow. □ 

Thus we have the following four characterizations of R(A, B ) which are equivalent to the origi¬ 
nal. 

(1) R(A, B) = sup {1/S (/) : /U = 0, /\ B = 1} 

(2) R(A,B) = 1/Sff) where f G £(V) is the the unique function with <f\ A = O,0| B = 1, 
A0(p) = 0 for p ^ A U B. 

(3) f?(A5) = inf {E(J) : flu x(A,B,J) = 1}. 

(4) R(A, B) = E(/) where / G £(E) is the unique energy dissipation minimizing unit flow 
from A to B. 


Proposition 5.2. Let G = (V, E) be a finite weighted graph with conductances c p q and suppose 
G = (V, E) is a weighted graph such that V = V and E C E with conductances c p , q = c pjj 
whenever (p, q ) G E. Then 

S G {f)>S d {f) for all fe£{V)=£{V), 

and hence 

R g (A, B) < R d (A, B ) for all A, B C V = V. 


Proof This follows because S G (f) is a sum over the edges in E and Sg(f) is a sum over a subset 
of the terms. □ 


Proposition 5.3. If G = ( V,E ) is a finite graph with conductances {c pP1 } (p q)C - E . Define G = 

(V, E) where V = Vj ~ where ~ is the equivalence relation defined by p ~ q if either p.q G /l, 
p, q £ B or p = q. That is, ~ identifies all the elements of A and all the elements and B. If G is 
given conductances 


c [po],[i?o] 


J2 c p« 


pe[p 0 ] ge[q 0 ] 


then 


R G (A, B) = R G (a,b) 


where a, b G V are the equivalence classes containing the elements of A and B respectively. 


Remark 2. The definition of c can be interpreted as coming from KirchhofFs law that conductances 
in parallel add. 

Proof If 0 G £{V) is taken to be the harmonic function such that f\ A = 0, 4 >\b = 1, and A G (j)(x) = 
Oforalla; f AUB. The function^ G £(V) with0([p]) = <f>(p) is well defined because f is constant 
on the fibers of namely A and B. Further, p{a) = 0, (p(b) = 1 and it is elementary to show that 

Aq4>(x) = 0 for x G V, x a, b. Thus R G (A, B ) = S G ((j)) = S G (4>) = R G (a, b). □ 
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